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Abstract 

Asymptotic independence of the components of random vectors is a concept used in many applica- 
tions. The standard criteria for checking asymptotic independence are given in terms of distribution 
functions (dfs). Dfs are rarely available in an explicit form, especially in the multivariate case. Often 
we are given the form of the density or, via the shape of the data clouds, one can obtain a good 
geometric image of the asymptotic shape of the level sets of the density. This paper establishes 
a simple sufficient condition for asymptotic independence for light-tailed densities in terms of this 
asymptotic shape. This condition extends Sibuya's classic result on asymptotic independence for 
Gaussian densities. 
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1 Introduction 

The purpose of the present paper is to provide simple sufRcient conditions that ensure asymptotic inde- 
pendence of the components of random vectors whose probability distribution is described by a density. 
Standard criteria for checking asymptotic independence are given in terms of distribution functions (dfs). 
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However, these are not always available in an explicit form in the multivariate case, and they give little 
insight in what large samples from a distribution will look like. Often we are given a density in analytic 
form. For light-tailed densities, the data clouds give a good geometric image of the asymptotic shape 
of the level sets. Hence it is of interest to have conditions for asymptotic independence in terms of the 
shape of the level sets of the underlying density, or in terms of a limiting shape for data clouds. 

For vector valued data it is standard practice to plot the bivariate sample clouds for all component 
pairs. In our final result it is the asymptotic behaviour of the shape of these bivariate sample clouds, 
as the size of the data set increases, that determines asymptotic independence of the coordinates for the 
underlying multivariate distribution. 

Unimodal densities whose level sets all have the same shape are called homothetic. The decay along 
any ray then is the same up to a scale constant depending on the direction, and hence the concept of 
light and heavy tails is well-defined. This remains true if we only assume that the level sets have the 
same shape asymptotically. Our primary focus here is on light-tailed densities, but for additional insight 
we also include some results for the heavy-tailed case. Throughout the paper, we assume continuity of 
dfs and of densities. 

Intuitively, for bivariate data, asymptotic independence means that large values in one coordinate 
are unlikely to be accompanied by large values in the other coordinate. Situations with a low chance 
of simultaneous extremes are often encountered in practice, for example in applications which involve 
modeling environmental data (e.g. [17]) or network traffic data (e.g. [19]); see [25] for further references. 
It is a well-known result, dating back to 1959 (see [26]), that the components of a vector with a Gaussian 
density are asymptotically independent whatever the correlation. Asymptotic independence also holds for 
light-tailed elliptical densities (see e.g. [15]). As a further generalization we shall show that a vector with a 
continuously differentiable homothetic light-tailed density whose level sets are convex has asymptotically 
independent components. 

Theorem 1.1. Let D be a bounded open convex set in M"^, containing the origin, with a boundary 
(i.e., at each boundary point there is a unique tangent plane). Let c„ > decrease to zero such that 
Cn+i/cn — > and let rn be positive reals such that r„+i/r„ — > 1. IfK is a random vector in R"^ with a 
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continuous probability density f whose level sets satisfy 

{f > Cn} = TnD n > no, (1.1) 

then the components ofH- are asymptotically independent. 

We shall prove more general results, Theorem 3.10 and 3.11 below. The generalizations we introduce 
are simple, but they result in theorems for which some extra terminology has to be developed. 

• The condition that the level sets all be of a given shape is replaced by the condition that the level 
sets can be scaled to converge to a limit shape D. Theorem 1.1 remains valid if we replace (1.1) by 
the limit relation {/ > c„}/r„ D. 

• Asymptotic independence holds if it holds for the bivariate marginals; in this spirit, our theorem 
imposes conditions on the projection D12 of the set D on the xi, a;2-coordinate plane. 

• Convex level sets are replaced by star-shaped level sets. 

• The condition of a smooth (C^) boundary is replaced by a condition which only affects the maxi- 
mum: the coordinatewise supremum of the points in D12 should not be a boundary point of Di2- 
Such a set is called blunt. 

• Our final result (Theorem 4.3) is on bivariate sample clouds. If these converge onto a blunt star- 
shaped set, asymptotic independence holds. 

The conditions in Theorem 1.1 above on the sequences (c„) and (r„) ensure that the density has light 
tails. If Cn+i/cn 1 then the density has heavy tails, and the same condition on the shape of the level 
sets implies asymptotic dependence. 

The paper is organized as follows. In Section 2 we discuss the concept of asymptotic independence. 
Section 3 introduces the class H of continuous multivariate homothetic densities with star-shaped level 
sets. The simple structure of these densities makes them a good starting point in our investigation of 
the relation between the shape of the level sets and asymptotic independence. Here we formulate our 
main results (Theorem 3.10 and 3.11) for light-tailed densities, and a counterpart (Theorem 3.12) for 
heavy-tailed densities. In Section 4 we introduce sample clouds, and state conditions for asymptotic 
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independence in terms of their limiting shape (Theorem 4.3). Section 5 provides various examples. The 
appendix clears up a number of minor points in the main text, by providing supplementary results and 
counterexamples. 

2 Asymptotic independence 

In this section we discuss asymptotic independence for multivariate distributions. We begin with a 
heuristic approach in the bivariate setting. This will clarify the significance of the concept for risk 
management. We discuss the commonly used criteria for asymptotic independence, the relation with 
multivariate extreme value theory, and describe possible forms of asymptotic dependence. 

2.1 Heuristics 

In finance one is interested in the future value of stocks, say in one year's time. Let X and Y denote the 
future value of two stocks. Suppose the distribution of the pair {X, Y) is given by a df F with continuous 
marginals Fi and F2. One is concerned about the risk that the stocks have a low value at this future date. 
Let Xp denote the p-quantile of the first stock, X, and yp the p-quantile of the second. The probability 
that the values of both stocks lie below the p-quantile for some given small value of p e (0, 1) is F{xp, yp). 
A risk averse investor would like this probability (of simultaneous loss) to be small compared to the 
probability p of a loss in either of the stocks. Future extreme low values of the two stocks are said to be 
asymptotically independent if F{xp,yp) = o{p) for p — > 0. Asymptotic dependence will increase risk for 
a portfolio containing these two stocks. 

A different example of risk is presented by the yearly maxima for high water levels at particular points 
on the coast of Holland, Great Britain and the U.S., say IJmuiden, Harwich and New Orleans. Consider 
a data set of 200 observations stretching back to the beginning of the nineteenth century. We assume 
that the data have been standardized to offset tidal effects. Now for each of these locations pick the five 
largest values. This yields three subsets of five elements in the set of 200 years. One would suppose that 
there is considerable overlap between the five years selected for IJmuiden and for Harwich, since in both 
cases the cause is the same, a North Western storm in the North Atlantic forcing water into the funnel 
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formed by the West coast of Holland and Belgium and the South coast of Great Britain, and opening into 
the Southern Atlantic via the Channel. High water levels in New Orleans have a different cause. So sea 
levels in New Orleans and in IJmuiden should be asymptotically independent, and one would expect the 
corresponding five point subsets to be disjoint with high probability. In general, given a positive integer 
k and a sequence of independent observations Zi, Z2, . . . from a bivariate df F, one could look in the 
sample cloud of the first n points Zj = {Xi, Yi) at the k largest observations for the coordinate Xi and for 
the coordinate Y^. This yields two subsets of the index set {1, . . . ,n}. Let Pn{k) denote the probability 
that the two subsets have a point in common. In the case of asymptotic independence one would expect 
that for fixed k, Pn{k) — > as n — > 00. This is indeed the CclSGj clS shown in [12] (Proposition 2). 

A related way to understand extremal dependence is by looking at the probability of a record, Pn{^), 
the probability that the coordinatewise maximum of a sample of n points is given by one of the sample 
points. It is shown in [13] (Theorem 2) that this probability vanishes as n ^ 00 exactly when the 
underlying vector has asymptotically independent components. For the sake of completeness we give an 
alternative proof of this basic result in the Appendix; see Proposition A.l. 

2.2 Sibuya's condition 

In his seminal paper [26] on multivariate extremes, Sibuya shows that the components of bivariate vectors 
with normal densities, whatever their correlation, are asymptotically independent. For a vector Z = 
{X, Y) with df F and continuous marginals Fi and F2 he introduces a function P by 

P{F,{x),F2{y)) =¥{X>x,Y>y} = 1 + F{x,y) - F,{x) - F^iy). (2.2) 

The function P is well-defined and continuous on the unit square. We can now give Sibuya's condition 
for asymptotic independence (see Theorem 2 in [26]). 

Definition Let Z = (X, Y) have df F with continuous marginals. The components X and Y are 
asymptotically independent if the function P introduced above satisfies 

P{1- s,l- s) = o{s) s>0,s^O. (2.3) 

By abuse of language we also say that the vector Z or the df F is asymptotically independent. 
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Independence of a bivariate vector is not affected by the marginal distributions. Similarly for asymp- 
totic independence. It is preserved under coordinatewise increasing transformations. Since we assume 
that the marginals Fi and F2 are continuous, there is a unique function C on the unit square such that 
F{x,y) = C{Fi{x),F2{y)). The function C is known as the copw/a of the df F. It is a df on the unit square 
with uniform marginals. Sibuya's condition is a condition on the copula C since P{u, v) = 1+C{u, v)—u—v 
for u,v & [0, 1]. If the condition holds for a vector {X, Y) with a continuous df F, it automatically holds 
for any vector whose df is continuous and has the same copula as F. 

In this paper we consider asymptotic independence for maxima. For minima one would define asymp- 
totic independence similarly, in terms of the copula, by: C(s, s) = o(s) for s J. 0. 

Proposition 2.1. Suppose there exist a,b> such that 

P{1- as,l-bs)/s ^0 slO. (2.4) 
Then asymptotic independence holds, and the relation (2.4) is valid for all a,b> 0. 

Proof We may assume a < 6 by symmetry and a = 1 replacing s by as in the denominator. By 
monotonicity Sibuya's condition holds. A similar argument gives (2.4) for any positive a,b. ^ 

For a vector {X,Y) with marginal dfs Fi and F2, Sibuya's condition may be formulated in terms of 
conditional quantile exceedances as 

Xu{X, Y) := limP{X > (g) | Y > F^{q)} = 0, (2.5) 

where F^{q) := infja; G M | F{x) > q} for q e (0,1) denotes the (minimal) g-quantile of F. The limit 
A(7 e [0, 1], if it exists, is known as the upper tail dependence coefficient. So X and Y are asymptotically 
independent if and only if their upper tail dependence coefficient is zero. 

Sibuya's condition is simple, but the formulation in terms of survival probabilities is inconvenient. 
The quantiles in (2.5) may be hard to determine since this amounts to computing the inverse of the dfs. 
There is a simple criterium in terms of sums: 

Proposition 2.2. Let {Xi,X2) have a df with continuous marginals Fx and F2, and let 1 — i^j(<j„) = 1/n. 
If Sn '■= nP{Xi + X2 > tin + t2n} ^ then Xi and X2 are asymptotically independent. 
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Proof Sibuya's function P satisfies nP{l — 1/n, 1 — 1/n) = nF{Xi > fi„,X2 > t2n} < ^ 0. This 
gives Sibuya's condition (2.3). ^ 

Below we give criteria in terms of continuous curves x{t) = {xi{t),X2{t)), t > 0, for which Fi{xi{t)) 
and F2{x2{t)) tend to one as t ^ oo. 

Proposition 2.3. Let {Xi,X2) have df F with continuous marginals Fi and F2. The components Xi 
and X2 are asymptotically independent if and only if for any e > there exists a continuous curve x(t), 
t > 0, such that Pi{t) = F{X > Xi(t)} is positive and vanishes for t ^ 00 for i — 1,2, and such that 

¥{Xi > Xi{t),X2 > X2{t)}/Pi{t) < e t>t„ i = 1,2. (2.6) 

Proof Assume asymptotic independence. One may choose Xi{t) and X2{t) continuous and increasing 
such that pi{t) = P2{t) for all t. Then (2.3) gives (2.6). Now assume (2.6). Let u G {0,Po] where 
Po = mm{pi{to),P2{to)}- By symmetry we may assume pi{t) = u< P2{t). Then 

P(l -u,l-u)< F{Xi > Xi{t),X2 > X2{t)} < epi{t) = eu. 

This holds for all u G {0,po]. So Sibuya's condition is satisfied. ^ 
Asymptotic independence is preserved under quite severe deformations of the distribution. 

Corollary 2.4. Let C be a convex open cone in M** and let f and g he probability densities which are 
positive on C, vanish off C and for which the quotients f/g and g/f are bounded on C. Let Z = 

{Zi, . . . , Zd) have density f and let X = {Xi, . . . , Xd) have density g. If Z\ and Z2 are asymptotically 
independent, then so are Xi and X2 ■ 

Proof The inequality / < Mg implies by integration that the same inequalities hold for the univariate and 
bivariate marginal densities and tail functions. So P{Xi > x\{t),X2 > X2{t)} < MF{Zi > x\{t),Z2 > 
X2{t)}. Similarly the inequality g < Mf gives an inequality with the constant M for the tail probabilities: 
V{Zi > Xi{t)} < MV{Xi > Xi{t)} for i = 1,2. Now use (2.6) for {Zi, Z2) to establish the same relation 
for {X\,X2). The extra factor has no effect in the limit. % 

The concept of asymptotic independence has been refined by looking at the rate at which P(l — s,l — s) 
vanishes for s ^ 0, or more generally by looking at the behaviour of the survival function P(l — u,l — v) 
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for u, u — > 0. See [25] or [23]. This second order theory is called hidden regular variation. We shall not 
treat this subject in our paper. Since our interest is in multivariate densities rather than dfs, we give 
conditions on the density which ensure asymptotic independence of the components. So the assumption 
of continuous dfs is not restrictive. We shall also assume that the densities are continuous on a convex 
cone and vanish outside this cone. In our context, this cone will typically be either the whole space R'^ 
or the open positive orthant (O,©©)''. 

2.3 Asymptotic independence and multivariate extreme value theory 

Asymptotic independence has to do with extremes, more precisely with bivariate maxima. Knowledge of 
multivariate extreme value theory is not indispensable for understanding asymptotic independence, but 
it will help to better understand asymptotic dependence. 

For multivariate extreme value theory we have to assume that each of the marginals of the multivariate 
df F lies in the domain of attraction of a univariate extreme value limit law (see e.g. Section 0.3 in [24] for 
a definition). As above let Zi, Z2, . . . with Z„ = (X„, F„) be independent observations from the bivariate 
df F, and write Z^" for the nth coordinatewise partial maximum. For simplicity assume the marginals 
are equal with the standard Frechet df Fi{t) = H{t) = e~^/* on (0, 00). The normalized marginal maxima 
Un = X^"/n and Vn = Y'^'^/n again have the Frechet df H by the scaling property H"[nt) = H(t). The 
scaled bivariate maximum W„ = jn has df Gn(w) = F"(nw). Suppose (j„ converges weakly to a 
limit distribution G, known as a multivariate extreme value distribution or a max-stable distribution. The 
limit vector W — {U, V) lives on (0, 00)^, and the components U and V have a Frechet law. Asymptotic 
independence for the df F is equivalent to independence of the components of the limit vector (see e.g. 
Proposition 5.27 in [24] or Theorem 6.2.3 in [8]). 

The convergence F"-{nw) G{w) becomes easier to handle if one takes logarithms. Write G(w) = 
g-/{(w) gjjjQg _ logF(nw) is asymptotic to 1 — F(nw) for w > we may write the limit relation as 

n(l - F(nw)) ^ i?(w) w € A" = [0, oo)^ \ {(0, 0)}. (2.7) 

The left hand side i?„(w) — n{l — F{nw)) is a df of the measure pn = n7r„, where 7r„ is the probability 
distribution of the vector Z/n. So pn is the mean measure of the sample cloud iV„ = {Zi/n, . . . , Z„/n} 
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and i?n(w) = pn{[0, w]'=). Pointwise convergence G„ ^ G on [0, oo)^ implies pointwise convergence -R„ 
RonX and vague convergence pn p- One can prove weak convergence pn p on A'\ [0, w] for any w e 
(0, oo)^. It follows that the scaled sample clouds Nn converge in distribution to a Poisson point process 
N on X weakly on the complement of centered disks. The Poisson point process N is of interest since it 
gives an asymptotic description of the large vectors in the sample cloud. Moreover, the coordinatewise 
maximum W„ of the scaled sample cloud converges in distribution to the coordinatewise maximum 
W of AT as the sample size goes to infinity. 

Sibuya's condition holds precisely if R{u, v) = 1/-^ + (Indeed n(l — Fi{nt)) — > 1/i for i = 1, 2, and 
the survival function F{x,y) = P{X > x,Y > y} = P{Fi{x), F2{y)) satisfies nF{nw) ^ by Sibuya's 
condition and Proposition 2.1.) In this case the measure p has density zero on the open quadrant (0, oo)^ 
by differentiation (see also Proposition 5.24 in [24]). So p is the sum of two measures, one on the 
positive horizontal and one the positive vertical axis, both with density The point process N thus 

is the superposition of two Poisson point processes on these halfaxes, and the two point processes are 
independent since the halfaxes are disjoint. This yields a simple description of the behaviour of large 
sample clouds: Asymptotically there is no relation between very large observations in the horizontal 
direction, and very large observations in the vertical direction. 

So far we have looked at asymptotic independence for bivariate distributions. Unlike independence, 
for a multivariate df asymptotic independence holds if it holds for the bivariate marginals; see e.g. 
Remark 6.2.5 in [8]. Here is the argument. For simplicity assume F has standard Prechet marginals, and 
the multivariate maxima converge in distribution. The scaled sample clouds Nn = {Zi/n, . . . , Z„/n} then 
converge in distribution to a Poisson point process N on X = [0, oo)'* \ {0} with mean measure p weakly 
on the complement of any centered ball. If Sibuya's condition holds for the bivariate marginals, then the 
projection pij of p on the Xi,Xj-p\a,ne will live on the two positive halfaxes in this plane. It follows that 
p lives on the set of points in [0, oo)'' \ {0} which have at most one positive coordinate, the union of the 
d positive halfaxes. The restrictions of N to these halfaxes are independent (since the d positive halfaxes 
are disjoint). Hence the multivariate extreme value limit vector has independent components. 
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2.4 Asymptotic dependence 

Asymptotic dependence is an ambiguous term. Logically it means the absence of asymptotic indepen- 
dence. We shall usually interpret it in a more constructive manner to mean the existence of a max-stable 
limit law G = e~^. This in turn imphes convergence in distribution of the normaUzed sample clouds 
to a Poisson point process N, whose mean measure p has di R = — log G, as sketched above. The tail 
dependence coefficient Xu in (2.5) gives only a very restricted view of extremal dependence. As pointed 
out in Section 8.2 in [21], a positive coefficient does not imply that the underlying distribution belongs to 
the maximum domain of attraction of some extreme value limit law. The exponent measure p provides 
a much more informative description of asymptotic dependence. Points of N which do not lie on one of 
the axes denote a simultaneous occurrence of very large values in two or more coordinates in the corre- 
sponding point of the normalized sample. Thus the point process N (or equivalently p) gives a complete 
description of how the extreme upper order statistics in the different coordinates are linked. 

3 Densities and level sets 

The aim of this section is to give conditions in terms of the density which will guarantee asymptotic 
independence. Two aspects of densities play an important role in our analysis: the shape of the level sets 
and the tail behaviour. We consider densities which are completely specified by just these two quantities 
- a shape for the level sets, and a decreasing function governing the rate of decay of the tails. In many 
cases of practical interest the shape is a bounded open convex set, containing the origin. The density 
then is continuous if and only if the decreasing function is. We shall also consider such unimodal densities 
which vanish outside the positive orthant. 

3.1 nomothetic densities 

nomothetic densities are densities whose level sets are scaled copies of a given open set D: 

{/ > c} := {x e M'' I /(x) > c} = rci* < c < co := sup/. (3.8) 



3 DENSITIES AND LEVEL SETS 11 

We assume that Z) c M** is a bounded open star-shaped set. A star-shaped set has the property that 
with any point x it contains all points rx for < r < 1. Assume the set D contains the origin. If each 
ray intersects the boundary dD in one point, then the set D can be represented using a gauge function 
n/j : R** — > [0,oo) that satisfies (i) n£)(fx) = tnoi^) for t > 0, x e K'' (homogeneity property), and (ii) 
D = {x gR'^ I nD(x) < 1}. The conditions on the set D and the continuity of / ensure that the gauge 
function is continuous. If the set D is convex, then so is the gauge function. If in addition the set is 
symmetric, —D = D, then the gauge function is a norm on W^, and the set D is the open unit ball in this 
norm. For any bounded open star-shaped set D the sets nD, n > 0, form an increasing family. Their 
union Z?oo is an open cone. It is a proper cone if the origin is a boundary point of D. One may then still 
define the gauge function no- This now is a function on the open cone £>oo- It is continuous on the cone 
if each ray in the cone intersects the boundary of £> in a single point, see Proposition 3.2 below. 

Gauge functions allow us to give analytic expressions for homothetic densities with the same ease 
with which one handles spherically symmetric densities. Continuous homothetic densities have the form 
/(x) = /o(n£)(x)) for a decreasing continuous function /o on (0,oo). This formula holds for all x for 
which the ray through x intersects the set D. 

Definition 1. Let Vd denote the class of all bounded open star-shaped sets D c M'', for which the 

cone Doo — IJ^^ nD is convex, and for which the gauge function no is continuous on this cone. A 
density f on M"^ belongs to the class 'H{D) if the shape D belongs to the class V^, and if f is of the 
form /(x) = /o(n£)(x)), where the density generator /o : [0, oo) — > [0, oo) is decreasing, positive and 
continuous. We set f = outside the cone D^o on which the gauge function is defined. 

It is apparent from the above definition that densities in H arc (star) unimodal since all the level sets 
{/ > c}, c G (0,co) are star-shaped; cf. Section 2.2 in [9]. Typical examples of densities in Ti. to keep in 
mind are the multivariate centered normal densities and, more generally, elliptically symmetric densities, 
discussed for example in [10]. See also Example A.l. In certain applications, elliptical symmetry may be 
too restrictive. Densities in Ti give the flexibility to model directional irregularities present in the data 
clouds, and to handle distributions on the positive orthant. The regularity conditions exclude pathological 
sets; see Example A. 2. 
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Before we proceed by looking at properties of these densities, let us review some related classes of 
models proposed in the literature. The £p-spherical densities ([22]) extend the class of spherical densities 
by allowing level sets to be balls in £p-norm for any p > 1. A further generalization is given by so-called 
v-spherical densities ([11]), where the scale function v plays the same role as the gauge function no 
defined above. In fact, our class W(£>) is a subclass of the t;-spherical densities in that we restrict level 
sets to be bounded and star-shaped, and / to be continuous on D^. In a recent paper [1], the authors 
advocate to study densities in terms of their contours. 

For a given shape D, what conditions does the density generator /o have to satisfy in order that the 
function /(•) = /o(nD(-)) is a probability density on M''? By regarding the set below the graph of / as a 
pile of thin Z)-shaped slices we obtain the following partial integration result: 



Observing that \sD\ = s'^\D\ and letting t tend to infinity, we obtain the condition (cf. Equation (5) in 



The class of unimodal densities introduced above is invariant under linear transformations. If the 
vector X has density / G Ti-iD) then the vector Y = AX. has density g e Ti{E), where E is the image of 
the star-shaped (convex) set D under the linear transformation A, and hence is also star-shaped (convex). 
A nice illustration of this invariance is the extension of spherical distributions to elliptical ones. 

For densities in TL the distinction between light and heavy tails is crucial for asymptotic independence. 
It is determined by the behaviour of the generator /o at infinity. Let us recall the definitions of regular 
and rapid variation. 

Definition 2. A measurable function h : (0, oo) — > (0, oo) is regularly varying at oo with exponent 6, if 
for a; > 



P{X &tD}= I fo{nD{^))d^ = fo{t)\tD\ + \D\ f s''\dfo{s)\ = f fo{s)d\sD\, 



t > 0. 



(3.9) 




[11]): l = d\D\J^s^-'fo{s)ds. 



lim h{tx)/h(t) 



6* e M; 



(3.10) 



= X 



ifO = 0, then h is called slowly varying; h is rapidly varying at oo if 



oo, < a: < 1 



lim h{tx)/h{t) = < 



(3.11) 







x>l. 
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If the generator of a density / = foino) varies rapidly, the density has light tails. In our terminology 
/ has light tails precisely if /o is continuous, positive, and strictly decreasing on [0,oo), and if there is a 
strictly increasing sequence r„ > such that rn-\-i ~ r„ and /o(7'n)//o(''ra+i) oo. 

Rapid variation of /o allows us to give strong inequalities for the measure /x with density / = /o('^^>)• 
Proposition 3.1. Let jjL have density f e 'H{D) with a rapidly varying density generator fo- 

(i) For any e > 

n{rD'') « ii{rD\e-^rD) r ^ oo. (3.12) 

(ii) For any non-empty open set U C D 

nirD") « n{rU) r ^ oo. (3.13) 

Proof Rapid variation implies that /o(e*t) < fo{t)/M eventually for t — > oo, and hence the rings 
i?„ = e^^+hn \ e"hD, n > 0, have measure n{Rn+2) < e'^'^^ iJ,{R„) / M . The rings i?„ are disjoint, 
and their union is the complement of tD. On summing the odd and the even terms we find with 
r] = e'^'^y{M -e'^'^^y. 

OO 

fi{e^HD') = J2HR2n) + li{R2n+i)) < viK^o) + f^{Ri)) = w{e''hD \ tD). 

n=l 

This gives (3.12). The integral over a thin ring is much larger than the integral over the set outside the 
ring. Formally, take 6 = e/2 and r = t — e to obtain (3.12). To prove (3.13) take a non-empty open subset 
Uo C U whose closure lies in D. Then e'^'^Uo C .D if e is small, and hence by rapid variation of /o as 

above the infimum of / over rUo is much larger than the suprcmum of / over the ring R = rD \ e~'^rD. 
Since |i?|/|rJ7o| is a constant we conclude that n{R) « ii{rUo), and (3.13) follows from (3.12). ^ 

We now give some extra details on star-shaped sets. 

Proposition 3.2. Let D be a bounded open star-shaped set. Suppose for each non-zero vector x G D 
there is one positive real ro such that rox lies on the boundary of D. Then the gauge function is 
continuous on the cone Doc = Un^-^- 

Proof Let B denote an open unit ball. If D contains the origin, it contains a ball eB, and by homogeneity 
ud < S on SeB. So the gauge function then is continuous in the origin. Now suppose it is discontinuous in 
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a point p outside the origin. We may assume that ud (p) = 1 and that there is a sequence p„ e -Doo such 
that p„ ^ p and nz5(p„) ^ c ^ 1. Let r G (0, 1) be close to 1. Then noirpn) rc, and rp G D impHes 
rpn & D eventually, hence rc < 1 and so c < 1 since r is arbitrary. This implies p„ G £> eventually, and 
this also holds for sp„ for 1 < s < 1/c. Hence sp is a boundary point. But so is p. This contradicts our 
assumption on the boundary of £>. ^ 

In general densities in ^{D) are not closed under projection, even if D is convex. If a random 
vector {Xi, . . . , X^) has density / e 'H{D), then the density of [X^, . . . , X^_i) need not be homothetic, 
and the univariate marginals need not even be unimodal; see Example A. 3. There are some exceptions. 
Projections of spherical densities are spherical, and if the level sets are balls in for some p G [1, oo] then 
this also holds for projections along the coordinate axes. The class itself is closed under projection. 
If -D is a bounded open convex set containing the origin, then so is the vertical projection E of D onto 
the horizontal hyperplane; if the origin is a boundary point of D it may be an interior point of E. 

Proposition 3.3. Suppose D e ©d- Let E be the vertical projection of D onto the horizontal hyperplane. 

Then E e V^-i- 

Proof Write z = (x, y) to distinguish the horizontal and vertical part of the vector z. It is clear that E 
is a bounded open star-shaped set. Moreover E^ is the projection of the cone £>oo, and hence an open 
convex cone. We have to prove continuity of the gauge function ue on E^o- If is not continuous, 
there exists a vector x and a sequence x„ — > x such that n_E(x) = 1 and nB(x„) < cq < 1 because 
E is open. Hence nD(x, y) > 1 for all y and there exist y„ such that nD{'x^,yn) < cq. Since D is 
bounded, the sequence (2/„) is bounded, and we may assume that it converges to some element yo- Then 
— > (x, 2/0)5 and continuity of njj implies n£)(x, t/o) < cq < 1. Contradiction. ^ 

3.2 Densities whose level sets are asymptotically star-shaped 

In this section we relax the condition that all level sets have the same shape D G T>d to the condition 
that the level sets, properly scaled, converge to a set £> e 2?^. We restrict attention to the light-tailed 
setting. 

Definition 3. Let D e V^. A positive probability density f on R'^ belongs to the set A{D) if there exist 



3 DENSITIES AND LEVEL SETS 



15 



sequences c„ > and r„ — > oo with c„+i /c„ — > and r„+i r„ s«c/i t/iaf /or any e > eventually 

e-^rnD C {/ > c„} C eV„£), n > no- (3.14) 

We write {/ > c„}/r„ D. A continuous positive function f is shape equivalent to f if its level sets 
satisfy (3.14). 

The sequences r„ and c„ determine a set of continuous decreasing functions r] which satisfy ?7(c„) r„. 
All these functions r/(c) vary slowly for c — > 0+ by the assumption that c„+i/c„ — > and r„+i/r„ — > 1. It 
is this set of slowly varying functions rather than the particular sequences c„ and r„ which are of interest. 

Proposition 3.4. If the slowly varying function rj above is strictly decreasing, defined on (0, c^] for some 
Cri > 0, and vanishes in c^, then the inverse function fo = rj^ is a continuous positive strictly decreasing 
function on [0, oo) which varies rapidly in oo, and f = foino) is shape equivalent to g, and for all e> 

e"'{/ >c}c{g>c}c e'{f > c} < c < c^. (3.15) 

Proof Rapid variation of the inverse function fo follows from Theorem 2.4. 7(i) in [6]. Shape equivalence 
holds since {/ > c} = r}{c)D and ?7(c„) ~ r„ by assumption. % 

Since r„ — > oo one may take a strictly increasing subsequence such that the asymptotic equality 
Tn+i ~ Tn remains valid. Take cq > Ci. Any continuous strictly decreasing function 77 on (0, cq] which 
vanishes in co and has the value r„ in c„ satisfies the conditions of the proposition above. So there are 
many continuous strictly decreasing functions /o on [0,oo) which vary rapidly in 00 such that foinn) is 
shape equivalent to g. 

For functions in A{D) the inequalities in Proposition 3.1 also hold. Moreover, they have the nice 
property that the d marginals gi oi g £ A{D) will lie in A{Di), where Z>, is the projection of D on the 
ith coordinate. 

Lemma 3.5. Let g £ A{D). Let pr for r > 1 denote the measure with density u g{ru)/g{rq) for a 
fixed non-zero vector q. Then for any open set U which intersects D 

Pr{D') « Pr{U). (3.16) 
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Proof Let / = foinn) be shape equivalent to g and define Hr to have density /r(u) = /(ru)//o(r). It 
suffices to prove (3.16) for the measures pr with density 5r(u) = g{ru)/ fo{r). Let V be a non-empty 
open subset of U which lies in the complement of SqD for some So > 0, and whose closure hes in D. For 
r > ro/6 we may apply the pointwise inequality (3.17) below and conclude that pr{U) > e~'^'^p,r{e'^V), 
where we choose e > so small that e^'^V C D. Then, using (3.13), iJ,r{e~^D'') « ij,r{e^V) and another 
application of the pointwise inequahty gives (3.16). % 

Proposition 3.6. Suppose X has density f G A{D). Let Y = ^(X) be a non-zero linear combination of 
the components of^. Then Y has a density g € A{J) where J is the open interval £,{D). 

Proof We may assume that ^ is the vertical coordinate. The condition /(e'x) < f{x)/M for ||x|| > ro 
implies by integration over horizontal hyperplanes g{e'^y) < e^'^~^'>^ g{y) / M for |y| > ro. This gives 
rapid variation. The thin tails of / ensure that g is continuous. Let £,{D) = {a,b). We claim that 
the average of g over r{a,a + ea) and over r{b — €b,b), say Oi and hi, is much larger than over both 
r(a — ea, a) and r{b, b + eb), say a2 and 62, as r — > 00. It suffices to show that P{ra <Y<ra + era} and 
¥{rb — reb < Y < rb} are much larger than P{X s rD"}. This follows from (3.16) since the horizontal 
slice {ra < Xd < ra + era} contains the open set U = rD n {xd < ra + era}, and P{X e U} is much 
larger than P{X G rW^}. A similar argument holds for the strip {rb — erb < xa < rb}. Let c lie between 
ai A 02 and 61 V 62- Then g > c on [ra — 2era, rb + 2er6] and g < c holds off [ra — 2era, rb + 2er6]. Hence 
e~^'^rJ C {g > c} C e^^rJ holds for small e for sufficiently large r. ^ 

It is not known whether a similar result holds for the bivariate marginals of /. 

One can also define shape equivalence for functions with heavy tails. In that case shape equivalence 
is the same as asymptotic equality in infinity, see Proposition A. 2. For light tails the behaviour of the 
quotient of two shape equivalent functions may be very erratic. Example A. 4 in the Appendix exhibits 
some functions which are shape equivalent to the bivariate Gaussian density. 

The level set inequalities in (3.15) imply the pointwise inequalities: 

/(e'x) < ff(x) < /(e-'x) X G r,D'. (3.17) 
This shows that a function g = qf is shape equivalent to f — foino) for g = e"^ and /o = e~'^ if for every 
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e > 

(p{r) - (fiie'-r) < ipi-x.) < ip{r) - ip{e-'-r) ic€drD, r>re. (3.18) 
Rapid variation implies that the left side goes to — oo and the right to oo. Hence we find: 
Proposition 3.7. g € -4(-D) is shape equivalent to g = qg if q is positive and logq is bounded. 

For Weibull-like functions of the form / = fo^no) where /o = e~'^ for a continuous strictly increasing 
function which varies regularly with positive exponent there is a simple alternative description of shape 
equivalence. 

Lemma 3.8. Suppose fo = e~'^ where f is continuous, strictly increasing and varies regularly in oo with 
exponent 6 > 0. Then fo varies rapidly in oo and g = e~^ e A{D) is shape equivalent to foino) if and 
only if^{'x.n) ~ (/5(nD(x„)) holds for all sequences ||x„|| oo. 

Proof Regular variation of (p gives (p{e'^r) — (p{r) ~ (e^^ — l)(p{r) and (p{r) — (p{e~'^r) ~ (1 — e^''^)(p{r) 
for r oo. The claim then follows from (3.18) with 'y = (f — t/j. % 

Proposition 3.9. Suppose g = is a continuous positive density on M''. Let there exist a function v 
on which is positive outside a hounded set and a non-zero vector q such that 

7(inU„)/7(t„q) tn^ 00, u„ ^ u, u G M*^. (3.19) 

Then there exists a set D & Vd containing the origin, a positive constant 6 such that u = n^jj, and a 
continuous strictly increasing function ip on [0, oo) which varies regularly with exponent 9 in oo such that 
7 is asymptotic to <p>{nD) in oo, and g is shape equivalent to foino) where fo = e~f . 

Proof Set a{t) = ^{tq). Then a varies regularly in oo with exponent ^ > 0. For any unit vector co the 
function t j{tu)) varies regularly in oo since j{tWn) / a{t) has a positive limit 1/(0)) for ci;„ — > a; and t — > 00. 
Uniform convergence on compact sets implies that ly is continuous on the unit sphere, and a{st)/a{t) s^ 
implies 'y{stu>)/a{t) u{w)s^. This proves that u{ru) = r^i'{u), and j{ru) a{r)i'{u) x a(r) for 
||ra;|| — > 00 since (3.19) implies uniform convergence on compact sets, and v is continuous and hence 
bounded on compact sets. Hence a(r) ^ 00 and (3.19) with u„ = gives j^(0) — and continuity of 
u implies 9 > Q. Set D = {v < 1}. Then = nf^. Let p G dD. Then t 7(tp) varies regularly 
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with exponent 6, and we may choose positive, continuous and strictly increasing and asymptotic to 
this function for t oo. Then 'y{rLu) ^ a{r)i'{uj) gives 7(w) ~ (p{nD{'w)) for ||w|| oo. Lemma 3.8 
then shows that g is shape equivalent to e~'^*^"°^ and since e"*^ varies rapidly if if varies regularly with 
exponent ^ > it follows that g G A{D). ^ 

Remark 1. Condition (3.19) is related to multivariate regular variation; see e.g. Section 5.4.2 in [24]. 

3.3 Criteria for asymptotic independence 

A random vector X with a spherically symmetric density /(x) = /o(||x||2) will have asymptotically 
independent components if the generator /o varies rapidly, and it will have asymptotically dependent 
components if /o varies regularly (see e.g. Theorem 4.3 in [15], or Proposition 3.2 in [14]). These results 
remain valid for / e A{B) where B denotes the open unit ball. One may replace B by certain bounded 
open star-shaped sets D, as will be shown in Theorem 3.11 and 3.12 below. 

With any bounded open set D one may associate the open intervals D, = {ai,bi), i = l,...,d, 
obtained by projecting D onto the ith coordinate. Then (a, b) is the smallest open box containing D, and 
b = {bi,. . . ,bd) = sup D is the coordinatewise supremum of all points in D (and a = (ai, . . . , a^) = inf D). 

Asymptotic independence depends on the bivariate marginal distributions. Hence we introduce the 
projections of D on the {xi, xj) coordinates. For l<i<j<dwe denote by Dij c the projection of 
D onto the two-dimensional space spanned by the unit base vectors and . The sets Dij lie in I?2 by 
Proposition 3.3, and Dij fits exactly into the rectangle {ai,bi) x {aj,bj). The cone generated by D^j is 
the projection of the cone D^o generated by D. 

Definition 4. The set D gT>2 is blunt if the point (6i, 62) = supD does not lie in the closure of D. 

A bounded open convex set D in K'' is smooth in the boundary point p if there is a unique hyperplane 
which contains p but which does not intersect D, the tangent plane to £> at p. For a planar set D this 
means that p is not a vertex. If the convex hull of D is smooth in all points then all bivariate projections 
Dij are blunt. 

We can now state our main results for asymptotic independence in terms of densities. 
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Theorem 3.10. //X has a light-tailed homothetic density f G 'H{D), and D is convex with a smooth 
boundary, then for any two distinct unit vector a and b the random variables a'^X and b'^X are asymp- 
totically independent. The result remains valid if the density of X is in A{D) . 

Proof First assume a and b are linearly independent. Introduce new coordinates such that a and b 
become the first two base vectors ei and e2. It suflaces to check that the assumption holds for vertical 
tangent planes, hyperplanes which project onto a line in the two-dimensional a;i , a;2-plane. The character- 
ization of D is geometrical and so it is preserved under linear transformations. The projection D12 in the 
new coordinates is also convex and smooth. Hence £'12 is blunt, and we may apply Theorem 3.11 below. 
If b = —a the bivariate distribution lies on the counterdiagonal, y = —x, and asymptotic independence 
is trivial by applying Proposition 2.3 with x{t) — {t, t) for i > 0. % 

Theorem 3.11. Suppose X has density g G A{D). If the bivariate projection D12 is blunt then Xi and 
X2 are asymptotically independent. 

Proof There is a simple analytic argument. Let supDi2 = {bi,b2). The sum Y = Xi + X2 has density 
g € A{J) by Proposition 3.6 where J has upper endpoint b = svLp{x + y \ {x, y) G -D12} and 6 < 61 +62 by 
bluntness. Now use Sibuya's condition for sums in Proposition 2.2. We give a more probabilistic proof 
in Section 4. ^ 

3.4 Criteria for asymptotic dependence 

We now give the counterpart to Theorem 3.11 for heavy-tailed densities. We are interested in the case 
where the partial maxima go to infinity in all coordinates. If the shape D lies in the negative orthant 
then the coordinatewise maxima converge to 0, if it lies in a negative coordinate halfspace {xi < 0} all 
partial maxima will lie in this halfspace too. We exclude these cases in the theorem below. 

Theorem 3.12. Suppose X has density f ~ foino) with D GVd and fo a continuous strictly decreasing 
positive function on [0, 00) which varies regularly in 00 with exponent — {X + d) for some A > 0. Assume 
for each coordinate i G {l,...,d} the set D contains a point whose ith component is positive. The 
components of X are asymptotically dependent unless D is contained in the set S of points with at most 
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one positive coordinate, a union of d+1 orthants. The partial maxima X^" may be scaled to converge 
in law to a vector W whose components have df V{Wi <t} — e"^"*/*)'^ for positive constants oi, . . . , a^- 
The exponent measure of W is the image under the map x i-^- x"*" = (xi V 0, . . . , V 0) of the excess 
measure p with intensity c/n^'^. One may take c = 1 by a suitable choice of the scaling constants for 
the maxima. 

Proof Choose a point wq in the cone on the boundary of D, and for r > 1 set hr{w) = /(rw) //(rwo) = 
foirnoi^))/ foil"). By regular variation for w e Doo, w 7^ 0: 

Convergence h^ ^ h holds uniformly on the intersection of the cone with any ring rB\eB, where 
B denotes the open unit ball. (The function I/ud is bounded on such sets.) By Potter's theorem 
(Theorem 1.5.6 in [6]) for any e > there exists such that 

fo{rs)/fo{r) < 2s'/s^+'^ r>r,,s>l. (3.20) 

This yields an integrable majorant for the convergence /i^ — > /i on \ B. Lebesgue's dominated 
convergence theorem implies that /ir —> /i in on Doo \ B, and because of the uniform convergence 
above also on D^ \ cB for any e > 0. 

Let p{r) be the finite measure with density hr, and choose r„ so that /9(r„) has mass n. Then /9(r„) 
is the mean measure of the scaled sample cloud Nn ~ {Xi/r„, . . . , X„/r„}, and pirn) p weakly on 
Doo \ cB implies iV„ => N weakly on Dqo \ e-B where N is the Poisson point process on Doo with intensity 
h. This tells us that the maxima converge. The measure p is an excess measure on M'' \ {0}: 

p{rA) = p{A)/r^ r > 0, A & Borel set in R*^ \ {0}. 

For halfspaces A = {xi > 1} or A = {xi < — 1} this relation also holds and implies that the marginals pi 
of p satisfy the same relation, and hence there exist non-negative constants cf such that 

Pi[r, 00) = c+/r^ Pi(-oo, -r] = c^/r^ r > 0. (3.21) 

For the coordinatewise maxima of heavy-tailed distributions it is convenient to work on the non-negative 
orthant, and replace the vector X by X+, where we use the continuous map {x\,. . . ,Xd) = x 1-^ x+ = 
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(a;i V 0, . . . , a;d V 0). We shall write and for the images of Nn and N under this map, and denote 
the mean measures by p+ and p+ . By assumption each component has a positive probability of being 
positive. Hence W„ := maxA^+ = max. Nn max A'' = maxA^+ =: W. Write i?(w) = oo, w]'^) = 
p+([0, w]'^), ^ w > 0. Then the limit distribution H of the coordinatewise maxima is 

P{W < w} = P{Ar([0, w]'^) = 0} = e-'''"([°'*l') = e-"(^). 

Using (3.21) the same argument gives the marginals Hi(t) = e~^''*'°°) = e~'^^/* with = . Since x"*" 
lies on a coordinate axis precisely ifxG S, the exponent measure p+ lives on the positive halfaxes if and 
only if £> C 5. 1 



4 Sample clouds 

In this section we look at the asymptotic behaviour of clouds of independent observations from a given 
light-tailed distribution as the number of data points in the sample approaches infinity. In particular we 
are interested in the limiting shape of these sample clouds under suitable scaling. Remark that sample 
clouds can be viewed as finite point processes with a fixed number of points. The motivation for looking 
at sample clouds is threefold. First of all, there is a relation between the asymptotic shape of the level 
sets of the underlying light-tailed density and the limit set onto which corresponding scaled sample clouds 
converge; see Proposition 4.4 below. Secondly, the point process approach will yield an intuitive proof of 
our main results. Finally, for sample clouds projection on the horizontal hyperplane is simple: just delete 
the last coordinate for each sample point. 

As before we consider a sequence of i.i.d. random vectors {Xi, X2, . . .} from a continuous distribution 

on with density /. Let Nn = {Xi/s„, . . . ,X„/s„} denote a scaled n-point random sample (or sample 
cloud) with the scaling constant s„ > 0, s„ — > 00 for n — > 00. Alternatively, for any Borel set A C K'' 

n 

Nn{A)=^lA{Xi/Sn). 
i=l 

The mean measure of Nn is given by p„ = n7r„, where tt,,, is the distribution of the scaled vector Xi/s„. 
The intensity of the n-point point process N^ is /in(u) = ns^/(s„u) for u e R'^. 
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Under a suitable choice of s„, the scaled observations 'X.i/sn from a density / e A{D) with high 
probability will fill out the closure of the shape set D in the sense of the following definition. 

Definition 5. Let E be a compact set in R'^ and fin finite measures. We say that the measures /U„ 
converge onto E if ^,i(p + cB) oo for any e-ball centered in a point p £ E, and if fj,n{U'') for 
all open sets U containing E. The finite point processes Nn converge onto E ifF{Nn{U^) > 0} for 
open sets U containing E, and ifF{Nn{p + cB) > m} — > 1, m > 1, e > 0, p € S. 

We call the set E in the definition above the limit set. In fact, the limit set, if it exists, is always 
star-shaped (see Proposition 4.1 in [16]). The following simple criterion is useful for checking convergence 
(in probability) of scaled sample clouds (see [4] for a proof). 

Proposition 4.1. If Nn is an n-point sample cloud from a probability distribution 7r„ on R"^, then Nn 
converges onto E if the mean measures fin = niTn converge onto E. 

The next theorem gives a sufficient condition for asymptotic independence of a distribution in terms 
of the limit set of the associated sample clouds. We first prove a lemma. 

Lemma 4.2. Let W„ = {Un, Vn) be the componentwise maximum of the sample Nn = {Z„i, . . . , Z„„} 
from the distribution Wn on R^. Suppose n7r„{(0, oo)^} — > 0, n7r„{u > 0} — > oo, and mrn{v > 0} ^ oo. 
Then the probability Pn = P{W e Nn} that the coordinatewise maximum is a sample point vanishes for 
n oo. 

Proof Since W„ lies in the positive quadrant or C/„ < or < 0, one finds 

P{W„ e Nn} < F{Nn{{0, oof) > 0} + ¥{Nn{{u > 0}) = 0} + F{Nn{{v > 0}) = 0}. 
These binomial probabilities all three vanish for n ^ oo. ^ 

Tlieorem 4.3. Let Xi,X2, . . . be independent observations from a continuous df F on M**. Let D be an 
open bounded star-shaped set which belongs to V^. Suppose there exist scaling constants Sn such that the 
scaled sample clouds Nn = {Xi/s„, . . . X„/s„} converge onto the closure of D. If the bivariate projections 
Dij of the set D are blunt then F is asymptotically independent. 
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Proof For the sake of simplicity assume d = 2. Let b = (61, 62) denote the coordinatewise supremum of 
D. Since D is blunt there exists a ^ > such that the shifted quadrant (e~^b, 00) and the set e^D are 
disjoint (see Figure la). Let n7r„ be the mean measure of A^„. Then n7r„(e~'^b, 00) < mTn{e^Dy 
and both n-Kn{D f\ {v > e~*6i}) and mrn{D n {u > e~*62}) go to infinity by Definition 5. The lemma 
above applied to the shifted clouds shows that the probability p„ of a record in the sample cloud Nn 
vanishes for n ^ 00. % 

In order to complete the proof of Theorem 3.11, we now establish a link between the asymptotic shape 
of the scaled level sets of a light-tailed density and the shape of the limit set of the associated sample 
clouds. 

Proposition 4.4. Let Xi,X2,... denote i.i.d. random vectors from a density g S A{D). Then the 
sequence of scaled sample clouds Nn = {Xi/s„, . . . ,X„/s„} converges onto the closure of the set D as 
n ^ 00 if the scaling constants Sn are chosen appropriately. 

Proof The density g is shape equivalent to / = foino)', see Proposition 3.4. The function /r(u) = 
/(ru)//o(r) also lies in W(-D), and by rapid variation and monotonicity of /o for any M > 1 and e > 
eventually fr > M on e~'^D and fr < 1/M off e'^D. By the pointwise inequality (3.17) the functions 
g'r(u) = g{rvL)/ foir) satisfy the same inequalities eventually if we replace e by 2e. The measure pr with 
density gr satisfies priD") « Pr{U), r — > 00, for any open set U which intersects D by Lemma 3.5. 
Choose Sn such that ns5^/o(s„) = 1 for n > tiq. Then ps^ is the mean measure of the sample cloud Nn, 
and from Proposition 4.1 the sample clouds Nn converge onto D since their mean measures do. ^ 



5 Examples 

This section illustrates the applicability of Theorems 3.10, 3.11 and 4.3 in determining whether a given 
distribution has asymptotically independent components. We also wish to see whether the conditions are 
sharp. 

For any open bounded convex set D in the plane whose closure contains the origin, the function e~"° 
is integrable and hence /(x) = coe""^'^) for suitable cq > is a probability density. Along rays it is an 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

(a) (b) 

Figure 1: Points at the edge of simulated sample clouds of size n = 10^ from (a) a bivariate normal 
distribution with mean zero and correlation p = 0.1 and (b) a bivariate meta-Cauchy distribution with 
standard normal marginals, in both cases restricted to (0, oo)^. Sample points are scaled by factor 
Sn = y/2 logn. The filled circles indicate the coordinatewise maxima for these samples. The boundaries 
of the corresponding limit sets E = {x^ — 2pxy+y'^ < 1 — for (a) and E = + + l > 3||(a;, 
for (b) are depicted by solid curves. The sample points inside the dashed curves are not displayed. 

exponential function. If U is uniformly distributed on D then it is simple to decide when the coordinates 
Ui and U2 are asymptotically independent, but asymptotic independence of the the coordinates of the 
vector X with density / is a different matter, even in the simple example where D is the intersection of 
a disk of radius r = 2 centered in (1, —1) and the open set above the diagonal. The light-tailed examples 
below are of a more general nature. 

Example 1. (Rotund-exponential densities) Let X have a continuous homothetic density / with 
convex shape D and generating function /q. If has a boundary and /o varies rapidly, the coordinates 
are asymptotically independent, and the sample clouds, properly scaled, converge onto the closure of the 
set D by Theorem 3.10 and Proposition 3.9. Now assume more: D is rotund, i.e. the boundary dD is 
with positive definite curvature in every point. Also assume that the generating function /o is asymptotic 
to a von Mises function e^"^ . Then / is a so-called rotund-exponential density, see [3], Sections 9 and 10. 
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If one zooms in onto a boundary point of D so as to distinguish individual sample points, the sample 
clouds converge to a Gauss-exponential point process: Nn => N vaguely on R**. The limit A'' is a Poisson 
point process with intensity g{u) = e~"'*e~'"i^ ^'^d-i)/2^(|27r)(''~^)/^ if one chooses the normalization 

appropriately. Weak convergence holds on all halfspaces {ud > co+ciUiH \-Cd-iUd-i}- The restriction 

of g to the upper halfspace {ua > 0} is a probability density. The corresponding vector has independent 
components. This vector is the limit of the high risk scenarios X^, properly normalized, where 
is the vector X conditioned to lie in the halfspace H, and H moves off to infinity in the sense that 
P{X e H} — > 0. These results remain valid under certain perturbations (if the density / = e"'^'"^^ is 
multiplied by a flat function L, see Section 11 in [3]). Such a perturbation does not affect the asymptotic 
behaviour of the exponent tp, but may affect the limit shape (take ip = log^(l + r) with = + 1/^ and 
L = with A = log(l + 2r + cx), c G [—1, 1]. If we choose c = c(r) = sin(loglogr) then L still is flat, 
but the shape of the level sets no longer converges). See Proposition 14.1 in [3]. So we see that under the 
extra conditions on the homothetic density / there are three alternative asymptotic descriptions of large 
sample clouds from this density. Convergence onto the closure of the set D describes the global behaviour 
of the sample clouds; weak convergence in the space X = [— oo, oo]*^ \ {— oo} to a Poisson point process 
whose mean measure is the exponent measure of an extreme value limit law (Gumbel with independent 
components); weak convergence to the Gauss-exponential Poisson point process A'' on certain halfspaces 
describes the local behaviour in boundary points of D. 

Example 2. (Skew- normal densities) A symmetric density g satisfies g{—x) = g{x). It may be 
transformed into an asymmetric density by multiplication with a positive continuous asymmetric function 
6 which satisfies ^(x) + 6'(— x) = 2. The skew-normal distributions SN{Q, a.) introduced in [2] have 
a density / which is the product of a centered Gaussian density with covariance f2 and the function 
X I— > 2$(a'^x), where $ is the standard normal df and a a non-zero linear functional. These densities 
are log-concave (since (— log$)" is positive) and hence have convex level sets. 

We claim that / € A{D) for a convex set D with boundary. Write / = e~'^. The function 
— log$(f) is asymptotic to for t — > — oo and vanishes for f — > oo. Let uq satisfy cx^uo > and 
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Figure 2: (a) Contours of the density of SN2{a; to) distribution with iv = 0.5 and a = (—1, 3); levels are 
indicated as powers of 10~^. (b) Asymptotic shape of the scaled level sets (bold curve); the dash-dotted 
straight line is given by equation a^u = 0. 

u^Q-^uo = 1. Then h{tuo)/t'^ 1/2 and 



h(tu) 
h{tuo) 



t oo. 



u^fi-^u for a^u > 

u'^fi-^u + (a^u)2 for a^u < 
It follows from Proposition 3.9 that there is a limit shape: D = {v < 1}, which is the covariance ellipsoid 
on the halfspace where oc^u is positive, and a flattened version of this ellipsoid on the complementary 
halfspace, see Figure 2. The set D is convex. To see that the boundary is C^, choose coordinates 
such that the underlying Gaussian density is standard, with spherical level sets, and then choose the 
vertical coordinate in the direction of a. In these coordinates D agrees with the unit ball B on the 
upper halfspace, and with the cylinder symmetric ellipsoid {x^ + • • • + + (1 + c^)x^ < 1} for some 
c > on the lower halfspace {xd < 0}. For boundary points p in the horizontal coordinate plane the 
tangent plane is vertical: p^x = 1. Theorem 3.10 applies. If X has a skew-normal distribution then 
the random variables ^i(X), . . . ,^rnPQ are asymptotically independent whenever the linear functional 
^1, ■ ■ ■ are linearly independent. Asymptotic independence of the skew-normal distribution has been 
partially proven in [18] using a direct analytic approach based on Sibuya's condition. 
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Example 3. Densities of the form /(x) = /o(||x||j,) for p > 1 have level sets which are balls in ip. Lower 
dimensional marginals have the same form but with a different generator, since the projection of the 
d-dimensional unit ball on the space spanned by the first m coordinates is the m-dimensional unit ball. 
The two-dimensional unit ball is blunt for all p G [1, oo), and hence vectors X with light-tailed densities 
/ as above have asymptotically independent coordinates. However for p = oo, the sup- norm, the unit 
ball is a cube C = (—1, 1)'', and the square is not blunt. The components of X arc still asymptotically 
independent, but we need extra work to prove this. The bivariate margins of /(x) = /o(||x||oo) have 
the same form with a different generator /o which still is continuous, strictly decreasing and rapidly 
varying in oo. (The cubic slices are replaced by square slices.) It suffices to consider bivariate densities 
f{x,y) = fo{\x\ V |y|) for continuous, strictly decreasing positive functions /o on [0,oo) which vary 
rapidly in oo. The marginals /j of /, i = 1,2, are equal by symmetry, and f2{y) = ^ufoiv) + 2-R(j/) 
where R{y) = /^°° fo{t)dt « yfo{y) by rapid variation of /q. Hence 2R{t) / f\{t) — > for t — > oo, and by 
I'Hospital's rule also the quotient P{X > t,Y > t}/T{X > t}. Thus Sibuya's condition holds. However, 
asymptotic independence need not hold if the level sets are only asymptotically cubic, see Example 2 in 



Another example showing that results for 'H{D) do not need to carry over to A{D) is given in the 
Appendix, Example A. 5. Our last example illustrates an asymptotically dependent distribution with a 
density in A{D) where bivariate projections of D are non-blunt. 

Example 4. Let Z have a bivariate t density with A > degrees of freedom. Transform the marginals 
to obtain a vector X with standard Gaussian components. The new density g is called a meta-t density; 
see [20] p. 193. It has normal marginals but the copula of the elliptic t distribution. The shape of the 
level sets of the density g converges to the symmetric subset D = {u\ + ^3 + A > (A + 2)|| (ui, U2)||L} of 
the square (—1, 1)^, see [4]. Figure lb shows a detail. The set D is not blunt, and the components of X 
are asymptotically dependent since those of Z are. 



[5]. 
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6 Conclusion 

We have explored conditions for asymptotic independence of the components of a multivariate random 
vector expressed in terms of the limiting shape of the level sets of the underlying density. A distinction 
had to be made between light and heavy tails. For light-tailed densities, the limiting shape of level sets is 
essential in determining whether asymptotic independence holds. In contrast, for heavy-tailed densities 
the (limiting) shape of level sets is irrelevant as long as the shape intersects the positive orthant. In the 
light-tailed case there is a simple sufficient condition for asymptotic independence of two components 
in terms of the corresponding bivariate projection of the shape. This subset of the plane has to be 
blunt. The more delicate question of the relation between shape and asymptotic independence when the 
bivariate projection is not blunt will be treated in a future publication. 

Asymptotic dependence is a basic concern in multivariate risk analysis. The light-tailed densities 
studied in this paper have the property that sample clouds will have the same shape as the level sets of 
the density asymptotically. For sample clouds persistence of the shape, as the number of sample points 
increases, opens the possibility of using the shape to construct densities over the whole space. This makes 
it possible to estimate probabilities of regions far out which contain only a few or no sample points. 
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A Appendix 

A.l Supplementary results 

Proposition A.l. Let Xi,X2,... be independent observations from the continuous df F on M^. The 
probability of a record amongst the first n observations goes to zero if and only if F is asymptotically 
independent. 
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Proof We may assume that F is a copula. Set c„ = nP{l — 1/n, 1 — 1/n) . Then the Poisson approximation 

gives a probability p„ = c„e~'^"e~(^~'^"^e~(^~'^"^ = Cnc'^""^ to the event: among the first n observations 
there is exactly one in the complement of [0, 1 — and that observation lies in the square (1 — 1/n, 1]^.. 
In case of asymptotic dependence Cfc„ — > c > for some subsequence, and hence the probability of 
a record in a sample of size fc„ will exceed ce'^~^/2 eventually. Conversely, asymptotic independence 
implies nP{l — M/n, 1 — M/n) — > for any M > 1, whereas the marginals satisfy nP{0, 1 — M/n) = M . 
As in Lemma 4.2 the probability of a record vanishes. % 

Proposition A. 2. Suppose f € 'H{D) and the generator fo varies regularly or satisfies O-variation 

r„+i ~ r„ ^ oo /o(r„+i) ~ /o(r„). (A.l) 
Then g is shape equivalent to f if and only if g ^ f. 

Proof Condition (A.l) implies that for every e > there exists ^ > such that /o(e^r) > fo{r) for all 
r > 0. Let n£)(x„) oo. Suppose g(x„) = c„ = fa{rn)- Then for any e > the point x„ eventually lies 
in the ring e^VnD \ e~^rnD on which / fluctuates by a factor at most e^. To show the converse, suppose 
fo satisfles (A.l) and let /(x„) = /o(r„) so that g{xn) ^ /(x„) = /o(r-„) ^ fo{rn+i) for r„+i ~ r„ ^ oo. 
Then for any e, ei > eventually 

(1 - ei)/o(eV„) < (1 - ei)/o(r„+i) < 5(x„) < (1 + ei)/o(r„+i) < (1 + ei)/o(e-V„), 

and since ei is arbitrary, we have /(e*x„) < g{xn) < /(e~^x„) for n > no as required for shape equivalence 
by (3.14). 1 

A. 2 Counterexamples 

This section contains counterexamples mentioned in the main text. 

Example A.l. A density / G 'H{D) may have spherical level sets without exhibiting spherical symmetry. 
This will be the case if D is an off-center ball with n£)(x) = ^ ||x||2 + /3^x — /3^x for some (3 gM.^. To 
be star-shaped the origin has to lie in D, or be a boundary point. In the latter case the set {/ > 0} is a 
ball (if {fo > 0} is a bounded interval) or an open halfspace. 
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Example A. 2. For any e > there exists a bounded open star-shaped set D which contains the origin, 
whose closure is the cube K = [—1, 1]'^, and such that the volume of D is small, \D\ < e. To construct 
such a star-shaped set take a dense sequence x„ on the boundary dK, and define U C dK as the union 
of open disks with center x„ and radius e„, where e„ ^ so fast that the area of ?7 is e/2. Now let D 
be the union of an open centered ball with volume e/2 and the set of all points ru with < r < 1 and 
uG i7. 

Example A. 3. Consider a continuous strictly positive density / on whose level sets {/ > c} = rcD 
all have the same shape. The set D is convex, even rotund, and the function c i—> rc is continuous and 
strictly increasing. The marginal densities are not necessarily all unimodal. 

Let /o be the uniform density on the tetrahedron Tq with two vertices in the horizontal plane, say ei 
and 62. The other two vertices arc e^/m and —63. Hero m is a positive integer to be chosen later. The 
marginal along the vertical axis along the base vector es has a continuous density go on (— 1, 1/m) which 
is parabolic on the interval (—1,0) and on (0, 1/m), vanishing in the endpoints of the interval and with a 
max;imum in the origin. The shifted tetrahedron T = To + e3/2 — (ei +e2)/8 contains the origin as interior 
point. Its vertical marginal density g is go shifted upwards over 1/2 and has its maximum in 1/2. The 
vertical marginal density g of the uniform distribution on the tetrahedron T/2 has its maximum in 1/4. 
The fair mixture of the uniform distribution on T and T/2 has a density / whose marginal g = {g + g)/2 
is not convex if m is large (since the left derivative of ^ in 1/4 is large). Now choose rotund sets £)„ 
converging to T and densities /„ converging to / which satisiy the conditions of the example. If infinitely 
many of the vertical marginals 5„ were unimodal then the limit g would be. We conclude that eventually 
gn is not unimodal. 

Example A. 4. In view of Lemma 3.8 we see that qg is shape equivalent to g for a continuous function 
q = precisely if |V'(x)| < xdl^ll) for a function x(r) << r^. Here are some examples of functions 
h = qg which are shape equivalent to the standard normal density g on the plane, for a continuous 
positive function q. One may take q to be one of the following functions 1 + |x| , (1 + r^)"* with m > 1, e'', 
gx-\v\^^^ where = x'^+y'^. These functions may be multiplied with a function like exp(sin7re^^ siuTre^") 
which fluctuates rapidly but is weakly asymptotic to a constant. The level sets {h > c} then will look 
locally like a shore with many small islands, and lakes, even though the sets are asymptotic to disks 
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{x^ + ]p < r^} with r = ^21og(l/c). 







Example A. 5. Let D be the open triangle with vertices (1,1), (—1,0) and (0,-1). It contains the 
origin. Let / € A{D) have convex level sets and be shape equivalent to g = e~"^ . It is possible that / 
is asymptotically independent. Suppose the function g has triangular level sets {g > e~*} = tD. Let / 
have level sets {/ > e~*} = Dt = tD\{x + y > 2t — Vi}. Then Dt/t agrees with the triangle D except 
that the extreme top has been sliced off. If we choose p{t) = {t,t) — {-Jt, \/f)/2, then 



where {h > e~*} = C + {t, t) for the halfspace C = {u + v < 0}. Let pt have density ht and let p have 
density h. Then pt ^ p weakly on [0,oo)^. Since p[0,oo)^ is finite and p gives infinite weight to the 
halfspaces {y > 0} and {x > 0}, Sibuya's condition holds by Proposition 2.3 with curve p{t), t> 1. 
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